Extra Dirac points in the energy spectrum for superlattices on single-layer graphene 
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We investigate the emergence of extra Dirac points in the electronic structure of a periodically 
spaced barrier system, i.e., a superlattice, on single-layer graphene, using a Dirac-type Hamiltonian. 
Using square barriers allows us to find analytic expressions for the occurrence and location of these 
new Dirac points in k space and for the renormalization of the electron velocity near them in the 
low-energy range. In the general case of unequal barrier and well widths the new Dirac points move 
away from the Fermi level and for given heights of the potential barriers there is a minimum and 
maximum barrier width outside of which the new Dirac points disappear. The effect of these extra 
Dirac points on the density of states and on the conductivity is investigated. 
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I. INTRODUCTION 

Graphene, a one-atom thick layer of carbon atoms, 
has been a topic of intense study since it's experimen- 
tal realization 1 in 2004. Interest in graphene results, 
in particular, from the prediction that carriers in it be- 
have as massless, chiral Dirac fermions, moving in a two- 
dimensional (2D) plane and described by a Dirac-type 
Hamiltonian. This model predicts unusual electronic 
properties such as the gapless electronic spectrum, the 
perfect transmission at normal incidence through any 
potential barrier, i.e., the Klein paradoxP^ which was 
recently addressed experiment ally^, the zitterbewegung 
recently verifiecP, etc., see Ref. [6] for a recent review. 

Motivated by all these properties, condensed matter 
physicists started extending the known properties of a 
two-dimensional electron gas (2DEG) in semiconductor 
materials to the relativistic 2D fermions (Dirac electrons) 
in graphene. Another particularly interesting system 
to consider is the application of a periodic potential to 
graphene, that is, a superlattice (SL), which under spe- 
cial conditions leads to collimation of electron beams^^l. 

As found recently in Ref. [TOl using a tight-binding 
formalism, the dispersion relation for such a SL can re- 
veal extra Dirac points at the Fermi leveffD. Close to 
the Fermi level the electronic properties of graphene are 
well described by the massless 2D Dirac equation. In 
two recent studies, Ref. [12] and [13j an exact condition 
was found for the emergence of extra Dirac points (zero 
modes) in the presence of a sinusoidal or square-wave SL 
potential. However, both studies are not able to describe 
the character and spatial distribution in k space of these 
new Dirac points, as they expand the spectrum for small 
k y . In this work we describe under which condition this 
is possible and also where these extra Dirac points arise 
in the electronic structure of massless Dirac fermions in 
single-layer graphene when a square-wave periodic po- 
tential is applied. Further, we analytically investigate 
the anisotropic renormalization of the group velocities 
at these new Dirac points, and find that the degree of 
the renormalization depends on the parameter u oc VqL, 



which is linear in the barrier height Vq and period L of 
the SL, in the sense that, v x > v y holds for u such that 
an extra Dirac point arises while for very high values u 
we have v x < v y . Moreover, we also consider the case 
of unequal barrier and well widths, not treated previ- 
ously, that results in a qualitatively different electronic 
spectrum. 

We organize the paper as follows. In Sec. ITT] we intro- 



duce our model. In Sec. Ill we investigate the emergence 
of the extra Dirac points, approximate the implicit dis- 
persion relation for small energies, and take a closer look 
at the group velocity near the extra Dirac points. Further 
we investigate the influence of the features of the spec- 
trum on the density of states and conductivity. Finally, 
we make a summary and concluding remarks in Sec. |IV| 



II. MODEL 
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Figure 1: Schematics of the superlattice potential V(x) con- 
sisting of square barriers. 

We describe the electronic structure of an infinitely 
large flat graphene flake by the nearest-neighbour, tight- 
binding model and consider solutions with energy and 
wave vector close to the K point. The relevant Hamilto- 
nian is % — vpa ■ p + IV , with 1 the 2x2 unit matrix. 
Explicitly % is given by 



u 



V —ivph(d x — idy) 

—ivph{d x + idy) V 



(1) 



where p is the momentum operator and vf ~ 10 6 m/s 
the Fermi velocity. In the presence of a one-dimensional 



2 



(ID) square-wave potential V(x), such as the one shown 
in Fig. 1, the equation (H — E)ip = admits solutions of 
the form ip(x)e tkyV with 



tp(x) 



ip(x) 



-i(f> 



-iXx 



(2) 



A = [(s — u(x)) 2 — k 2 ] 1 / 2 , tancj) = k y /X, s = sign(e — 
u(x)), e = E/vpfi, and u(x) = V{x)/vph] the parameters 
e and u{x) are in units of inverse length. 



A. Square-barrier superlattice 

We consider an infinite number of periodically spaced 
barriers, as shown in Fig. [I] with unit cell length L 
and barrier (well) width Wj, (W w ). It is convenient to 



introduce the dimensionless variables e — >■ eL, k v 



kyL, 



k x — > k x L, u — > uqL — VqL/vfH, x —> x/L, Wb —> Wb/L, 
and W w —> 1 — Wb/L. The wave function of this periodic 
system is a Bloch function and the transfer matrix T 
pertinent to it leads to an expression for the dispersion 
relation, see Appendix |A| For \k y \ < \e w \ and e w as in 
Eq. ^ , the transfer matrix T can be written 

asP 



T 



(3) 



(4) 



(5) 



then the dispersion relation becomes 
cos{k x ) = 5R{e- iA w}, 

with w given by 

w = e lXWb [cos KW b - iG sin AW b ] 
and 

e w = e + uWb, e h = e- uW w , G = (e w e b - k*)/\A, 
\=[el-kl}^, A = [4 (6) 

Writing the rhs of Eq. Q explicitly gives 

cos k x = cos \W W cos AWb — G sin XW W sin AWb- (7) 

Although this derivation is only correct for \k y \ < \s w \, 
Eq. ([7]) is also valid beyond this limitation, see Ref. US 
From Eq. ^ it can be seen that the dispersion relation 
possesses the symmetry property e — > — e for Wb -H- W w . 
The asymmetric spectrum is not unexpected because the 
symmetry of the potential about the Fermi level is lost 
for W b f 1/2. For W b = 1/2 we have 

7 A A A . A 

cos k x = cos — cos — — G sin — sin — , (8) 

where e w = e + u/2 and Eb = e — u/2. For this interesting 
case, the potential possesses particle-hole symmetry and 
the extra Dirac points originate at the Fermi level; we 
will show their arrangement, in k space, in Sec. Ill A 
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Figure 2: (Color online) (a) Valence and conduction bands of 
the spectrum of a SL with square barriers of width Wb = 1/2 
and height u = IQtt. (b) Contourplot of the conduction band. 



In Fig. [2ja) the spectrum resulting from Eq. ^ for 
equal barrier and well widths, i.e., for W w = Wb = 1/2, 
is plotted for u — IO71". As can be seen, the spectrum 
is symmetric about the Fermi level; there are two ex- 
tra Dirac points on both sides of the main Dirac point, 
and their velocities are renormalized. The anisotropic 
behaviour of the new Dirac cones is clearer in the projec- 
tion of the conduction band shown in Fig. ||b). Further 
details about the renormalization of the velocities will be 
given in Sec. |IIIC| 




Figure 3: (Color online) (a) Valence and conduction bands of 
the spectrum of a SL with barriers of width Wb = 0.4 and 
height u = 67r. (b) and (c): projection contours of the con- 
duction and valence band, respectively, on the (k x , k y ) plane. 



For unequal barrier and well widths the spectrum is 
shown in Fig. ^a) for W b — 1 — W w — 0.4 and u = 
The spectrum is no longer symmetric about the Fermi 
level, the two extra Dirac points are shifted in energy 
relative to the main point and their velocities are renor- 
malized. The location of the extra Dirac points will be 



investigated in Sec. Ill 
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III. ELECTRONIC STRUCTURE 



A. Appearance of extra Dirac points 

In order to find the location of the Dirac points we 
assume k x = 0, e = 0, and Wi — W w = 1/2 in Eq. Q. 
Then Eq. (8) becomes 

1 = cos 2 A/2 + [(u 2 /4 + fc 2 )/(u 2 /4 - fc 2 )] sin 2 A/2, (9) 



which has solutions for it 2 /4 — fc 2 = w 2 /4 + fc 2 or 
sin 2 A/2 = 0. For the first possibility k y = is the only 
solution and corresponds to the usual Dirac point. The 
second possibility leads to A/2 = jir with j ^ 0, because 
A = makes the denominator u 2 /4 — fc 2 = A 2 vanish and 
does not lead to a solution. For A/2 = jir we have 



A: 



y,j± 



= ± 



- 4j 2 7T 2 = ± 



2Hvf 



2J7TN 2 



(10) 



where we reinserted the dimensions after the second 
equality sign. As such, Eq. ( 10 ) describes the spatial 



arrangement of the extra Dirac points along the k y axis. 
Also, it clearly shows how many points we have at partic- 
ular values of it, namely 2 x (u mod 4ir), and where they 
are located in k space. Each time u becomes a multiple 
of 47T a new pair of Dirac points is generated for k y = 0. 
The condition j ^ gives us a threshold value of u — Air 
for the emergence of the first pair. The integer j denotes 
the jth extra Dirac point, so the outer extra Dirac points 
have j — 1 as they are generated first. 




Figure 4: (Color online) Slices of the SL spectrum along k y 
with k x = and u — 6n. (a) The solid red, dash-dotted 
green, and dashed blue curves correspond to barrier widths 
Wb = 0.5, 0.6, and 0.7 respectively. The thin black line is 
the curve on which the extra Dirac point, to the left of the 
main one at k y = 0, is located for various Wb- Only the 
new points to the left of the main one are shown, (b) As in 
(a) for fixed Wb = 0.5. The solid red, dot-dot-dashed black, 
dashed green, and dash-dotted blue curves are for different 
values of the barrier height such that u/2 = 2n, 3tv, 47T, and 
6tt, respectively. 

In Fig. Qa) we show slices of the SL spectrum along 
k y for k x = and u = 67r. The solid red, dash-dotted 
green, and dashed blue dashed curves correspond to bar- 
rier widths Wb = 0.5,0.6, and 0.7 respectively. The thin 
black line is the curve on which the extra Dirac points, on 
the left of the main one at k y = 0, are located for various 



Wb- In Fig. [4jb) we show slices of the spectrum along k y 
for k x = 0. The solid red, dot-dot-dashed black, dashed 
green, and dash-dotted blue curves are for different val- 
ues of the barrier height such that u/2 = 2ir, 3tt, 4tt, and 
67T respectively. For values of u/2 which are multiples of 
2-7T, new Dirac points are generated. Interestingly, if new 
extra points are to arise, the dispersion becomes almost 
flat along the k y axis at the Dirac point, i.e, collimation 
occurs. We will come back to this issue in Sec. IIII CI 

Unequal well and barrier widths. We return to the more 
general case of unequal well and barrier widths for which 
Wb 7^ 1/2. It is more difficult to locate the extra Dirac 
points which no longer occur at the Fermi level as seen 
from the green and blue curves in Fig. |3Ja) showing 
slices of spectra from Eq. Q for k x = 0. By means 
of the symmetry e —> — e for Wb ■<-> W w , we know the 
complementary plots for Wb — ^ 1 — Wb- As can be seen, 
the extra Dirac points shift mainly down (up) in energy 
as Wb increases (decreases). To find their coordinates 
(e, k x = 0,k y ) we assume sin(\W w )=sm(AWb) — and 
cos(XW w ) = cos(AWb) = ±1^. This gives (Appendix[B| 



tt 2 (f 



2^- 2W ^2u\Wl 



(j + 2mf 



W 2 



(e i)TO + uW h f - (j*/W w ) s 



1/2 



(11) 



where j and m are integers. This method also shows 
higher and lower crossing points if m =/= 0. In Fig. Qa) 
the extra Dirac points on the left, obtained with this 
method, are indicated by open circles and the thin black 
curve shows their trajectory in (E, k y ) space as the width 
Wb varies. For a particular u there is a minimal width Wb 
(and a corresponding maximal width Wb — > 1— Wi) below 
(above) which the various extra Dirac points disappear. 
In Fig. |4|a) the "Dirac cones" at these crossing points 
for m = are not only reshaped with a renormalized 
anisotropic velocity but, as shown by the blue dashed 
curve, the "extra Dirac point" is not at a local minimum 
(maximum) of the conduction (valence) band. 



B. Analytical expression for the spectrum for small 
energies e 

As the purpose is to have a closer look at the behaviour 
of the extra Dirac points and we cannot prohibit k y from 
being large, we expand Eq. Q for small energies, up 
to second order in e, and obtain the following explicit 
dispersion relation 



e±=± 



4|a 2 | 2 [fc 2 sin 2 (a/2) + a 2 sm 2 (k x /2)] 
fc4 a sina + a 2 u 4 /16- 2fc2 u 2 sin 2 (a/2) 



-l 1/2 



(12) 



witha = [tt 2 /4-fc 2 ] 1 / 2 . 

If we only need the behaviour of the spectrum near 
the K point (for small k x and ky), it suffices to make an 
expansion for small e and k y in Eq. (8), up to third order 



4 



in products of e and k y since this is the first order with 
an energy dependence. The result is 



2 cos k x — 2 -\- e — fc 2 sin 2 (w/4)/(u/4) 2 = 0. 
Then we solve for the energy e and obtain 

e ps ±[4sin 2 k x /2 + k\ sin 2 (u/4)/(u/4) 2 ] 1/2 . 



(13) 



(14) 



In Fig. [5]ja) we show e from Eq. ( 12 ) and compare it 
with the exact dispersion relation, for k x — 0, following 
from Eq. (Tsl) . The expansion (12) is rather good for low 



energies near the extra Dirac points; accordingly, we will 
use Eq. (12 1 to further assess their behaviour. 




Figure 5: (Color online) (a) Plot of the dispersion relation 
resulting from Eq. ( 12 1 (dashed blue curve) versus the exact 
one obtained from Eq. (ISJ) (solid red curve) for k x — and 
u — 6n. (b) The solid blue and red curves show, respectively, 
the zoomed-in plots of the solid red and dashed green curves 
of Fig. |4|b), i.e., for u = 4-k and u = 87r, respectively, and 
Wb = W w — 1/2. The approximation of these curves by Eq. 
(17 1 are the dashed curves. 



C. Anisotropic velocity renormalization 
at the (extra) Dirac point (s). 

The spectrum in the low-energy range consists of two 
kinds of valleys, one near the main Dirac point and the 
other near the extra Dirac points. Near the original Dirac 
point the spectrum is almost linear, perpendicular to the 
barriers, and zero parallel to them, whereas near the ex- 
tra Dirac points the situation can be reversed depending 
on the height of the barriers. 

Group velocity along the x-axis at the Dirac point. To 
compare with the collimation found by Park et a/P, we 
notice that in Fig. Qb), for the solid red, dashed green, 
and dash-dotted blue curves, corresponding to barrier 
heights which are multiples of 47r, the dispersion becomes 



more flat for small k y . From Eq. (14) we could already 



expect that, to order k 2 the k y dependence disappears 



for these values of u. Further, if we expand Eq. (12) in 



powers of k y we obtain 



e = sin(u/4)/(w/4) k y - (2/m 5 ) [m 3 cos (u/4) 
+4m 2 sin(u/4) - 128 sin 3 (u/4)] k y + 0{k y ), 



(15) 



which is linear in k y , for small k y , and the velocity be- 
comes 



v y /vF = de/dk 



v 



i(u/4)/(u/4). 



(16) 



In Fig. [6] the velocities of the Dirac point, in the x and y 
directions, are given by the j = curves. For u/2 = 2jir 
we have 



±kt 



, -2 2 
]j 7T 



0{kl) 



(17) 



which is cubic in k y for small k y 



If j and consequently u 
become larger the dispersion gets flatter. In Fig. [5jb) we 
plot Eq. (17) for u = 4-7T and u = 8n as dashed curves, 



which correspond, respectively, to the zoomed-in plots of 
the solid red and dashed green curves of Fig. |4|b), shown 
here as blue and red curves. 

Group velocity along the y-axis at the extra Dirac 



points. The dispersion relation (12) for the k y val 



ues of the extra Dirac points, determined by k y ^ ± = 

±[u 2 /4— (2J7T) 2 ] 1 / 2 , gives us an idea of how dispersionless 
the spectrum near these points is along the ^-direction. 



If k y - ± exists, Eq. (|12|) becomes 



327rVsm(|fc a |/2)/u 2 , 



(18) 



and the partial derivative^ of e with respect to k x is 

v x/ v f — de/dk x rj sign(k x )16n 2 j 2 cos(k x / 2)/ u 2 . (19) 

This means that for smaller j (the most distant extra 
Dirac points) the group velocity along the x-direction is 
strongly suppressed. Further, as u > Ajir must hold in 
order for k Vt j± to be real, \v x \ is smaller than 1 (= vp) at 
k x = 0. Only for the special values u = ijir, for which 
new Dirac points appear, we have \vj x \ = 1. 

Meanwhile the dispersion in the k y direction is also 
of interest. First, 



let us take k x — and expand the 
dispersion relation ( 12 ) for k y — k y ^ ± -C 1. To first order 
in this difference we obtain 

sn±[4k y l ± /u 2 ](k y -k y . t± ). (20) 

This gives the velocity v y at the extra Dirac points 

^ = de/dk y a 4fc y 2 ± /u 2 = 4[u 2 /4~4fir 2 ]/u 2 . (21) 

Since the coordinates of the extra Dirac points should be 
real, k y 2 ± is positive and smaller than u 2 /4 and we have 
v y < 1 (the outer Dirac points, for j = 1, show the largest 
Vy). This entails that both v x and v y are renormalized 
at the new Dirac points. A plot of the velocities of the 
extra Dirac points, in the x and y directions, given by 
Eqs (19) and (21), is shown in Fig. [6] As seen, for the 



extra Dirac points, Vj x , shown by the dashed red curves, 
starts from vp and decreases to zero with increasing u 
while Vj y (dash-dotted blue curves) , starts from zero and 
approaches vp slowly for large u. The thin solid black 
curves show the velocities for j = 1, obtained numerically 
from Eq. ([8]). As can be seen, the two curves match 
almost perfectly. 
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Figure 6: (Color online) Velocities vq x and vo y (dash-dot- 
dotted red and solid blue curves, respectively), vs u at the 
original Dirac point and Vj x and Vj y (dashed red and dash- 



dotted blue curves, respectively), given by Eqs ( 19 1 and (21 1, 
vs u at the extra Dirac points j = 1,2,3. The thin black 
curves for the j = 1 Dirac point, are obtained numerically 
from Eq. @. 



D. Density of states and conductivity 

Density of states. At zero temperature the density of 
states (DOS) D(E) is given by 




Figure 7: (Color online) The DOS, for u = 6tt, and W b = 0.5 
and 0.4 is shown by, respectively, the solid red and dashed blue 
curve. Stars and crosses placed near the peaks of the solid 
red curve (for negative energies) are ascribed, respectively, to 
saddle points in the spectrum, for k x = 0, and to minima for 
k x = 7r. The DOS without a SL potential is shown by the 
dash-dotted black curve. 



pendence is given by (3 = Hvp /ksTL = 20 (in dimen- 
sionless units). The solid red and dashed blue curves 
correspond to Wb = 0.5 and Wb = 0.4, respectively. 
The dash-dotted black curve shows the conductivity at 
zero temperature and in the absence of a SL poten- 
tial, a xx — <j yy = SfUq/A-k, with Bp — EpL/hvp and 
(To = e 2 /h. Notice that a xx is an oscillating function 



D(E) = J2^( E - E nk), (22) 

n,k 

with E the energy. We show the DOS in Fig. [7J for 
Wb = 0.5 (solid red curve) and Wb = 0.4 (dashed blue 
curve), as well as the DOS for graphene without any SL 
potential (dash-dotted black curve); the latter is given 
by D(e) = eDq/2it, with Do = L/hvp the amount of 
states per unit area and L is the period of the SL. The 
DOS shows an oscillating behaviour. The dips in it are 
located at the crossing points in the energy bands for 
k y = (e = nir), while the peaks marked by a star are 
ascribed to the saddle points between the crossing points 
for k x = and to the minima of the energy bands at the 
edge of the Brillouin zone, k x = ±7r, marked by a cross. 
For W b = 0.4 ^ W w the DOS (dashed blue curve) does 
not vanish at e = nor is it symmetric about this energy. 

Conductivity tensor a. The diffusive dc conductivity 
er„„ for the SL system can be readily calculated from the 
spectrum if we assume a nearly constant relaxation time 
r(Ep) = Tp. It is given bjfil] 

o-fj, u {Ep) = 6 — — v nll v nu f nk {\ - / nk ), (23) 

n,k 

with A the area of the system, n the energy band index, 
/i, v = x,y, and /„ k = l/[exp(P(E F - E nk )) + 1] the 
equilibrium Fermi-Dirac distribution function; Ep is the 
Fermi energy and /3 = 1/ksT. 

In Figs. [8]ja) and[8jb) we show, respectively, a xx and 
G yy for a SL with u — 67r, and the temperature de- 




Figure 8: (Color online) (a) and (b) show the conductivities 
a xx and a yy vs Fermi energy for a SL with u — 6n. The solid 
red and dashed blue curves are, respectively, for Wb = 0.5 
and Wb = 0.4 (W w = 1 — Wb). The dash-dotted black curves 
show the conductivities in the absence of the SL potential, 
<y xx = <y yy ~ EF<?o/4:ir. The inset in (a) shows the zoomed- 
out conductivity a xx , for Wb = 0.5, and the dashed lines are 
the conductivities of graphene in the absence of a SL but with 
a constant non-zero potential applied, — Vo/2 and +Vo/2, so 
that Ef is displaced by Vb/2 and u/2 — 3n. The inset in (b) 
is a zoom on a yy for small energies, (c) and (d) show a xx and 
<J yy , for Wb = W w =0.5 for different potential heights, such 
that u = 4-7T, 4.5-7T, 6-7T, 7.5-7T, and small energies. 

of the Fermi level and recovers a quasi linear behaviour 
similar to that in graphene without a SL potential when 



G 



ef is well above the barrier height, as shown in the inset, 
i.e. with the x axis displaced over the potential height or 
well depth, i.e., u/2 — 3tt. On the average a yy increases 
with e f and approaches the result without a SL for large 
energies. 

The oscillations in both a xx and a yy result from the 
motion of the Fermi level through the different SL mini- 
bands. Notice that for Wb — 0.4 the conductivities are 
asymmetric with respect to electron and hole conduction. 
In both cases, Wb = W w = 0.5 and Wb = 0.4, a xx shows 
dips at ep — nir, where energy band crossings occur in 
the spectrum for k y = 0. In the former case the DOS 
has dips occurring at the same energy values that are 
dominated by the same crossings. In the latter case we 
see that, unlike the DOS, a xx is almost unaffected by the 
extra Dirac points for low energies since the spectrum is 
almost flat near these points. Similarly, for Wb = 0.4 we 
see that the minimum in a yy is located at Ep ~ 1, that 
is, the energy value for which the two extra Dirac points 
occur in the spectrum. 

In Figs. [8^c) and [8^d) we show, respectively, a xx and 
a yy , for a SL with Wb — W w = 0.5, for different po- 
tential heights, such that u = 0, 47r, 4.5-7T, 67r, 7.57T, and 
(3 = hvF/ksTL = 20. Notice that the conductivity a yy , 
in the low-energy range (ef < 1), is lower than that in 
the absence of a SL potential while its slope increases 
as the potential barriers become higher. This is due to 
the extra Dirac points, that appear for larger potential 
heights, near which the velocity is larger along the y-axis. 
Notice that for Ep < 1 we have a xx > a yy as a result of 
the inequality v x > v y near the Dirac point. 



the Fermi energy and found an explicit expression for the 
k space behaviour of the extra Dirac points at the Fermi 
level. Using this expression we showed how the group 
velocities at the various extra Dirac points are renormal- 
ized in the x and y directions. We also quantified how 
dispersionless the spectrum is in the neighbourhood of 
a Dirac point along the y-direction and the emergence 
of new points at which the conduction and the valence 
bands touch each other. 

Finally, we obtained numerically the density of states 
(DOS), which exhibits an interesting oscillatory be- 
haviour and is reflected in the conductivity of the sys- 
tem. We found that the dips in the DOS, for symmetric 
SLs, are located at the touching points in the spectrum 
for k = 0, i.e., for s = nn. For asymmetric SLs these 
dips persist but extra dips due to the extra Dirac points 
arise. The conductivity <j xx was found to have dips at the 
same values for ef as the DOS while the main features of 
a yv in the low-energy range are due to the spectrum near 
the extra Dirac points. We notice in passing that the in- 
fluence of velocity renormalization on transport was not 
studied in Refs. HH and [DS] nor the modification of the 
extra Dirac points for unequal well and barrier widths. 
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IV. CONCLUSIONS 

We investigated the appearance of zero modes, touch- 
ing points at the Fermi level (extra Dirac points) in the 
spectrum of single-layer graphene in the presence of a ID 
superlattice (SL). The system was described by a Dirac- 
type Hamiltonian, and the SL barriers were square. 

In the general case of unequal well and barrier widths, 
there is no particle-hole symmetry and the extra Dirac 
points are no longer located at the Fermi level. We ob- 
tained an analytical expression for the position of the 
crossing points in the spectrum. The extra "Dirac cones" 
that appear at the various crossing points are reshaped, 
i.e., they are no longer circular symmetric and the slope 
is renormalized. For fixed height of the barriers, we found 
lower and upper bounds for the barrier and well widths 
for the occurrence of these extra Dirac cones. 

For a SL with equal well and barrier widths we com- 
plemented the investigations of Refs. H2] and Q21 which 
numerically demonstrated the emergence of extra Dirac 
points (zero modes). In doing so we found a simple an- 
alytical expression for the spatial distribution of these 
points in k space as well as a threshold value of the po- 
tential strength for their appearance. Further, we ap- 
proximated the dispersion relation for energies close to 



Appendix A: Dispersion relation for periodic 
systems 

The wave functions in the regions before and after the 
barrier, labeled, respectively, by j = 1 and j — 2, can be 
written as tpj(x) = £lj(x)Aj, with 



n = 



A = 



(Al) 



Since the wave function of the entire periodic system is 
a Bloch function and the transfer matrix T connects the 
regions before and after the barriers, we have 



^(0), A 1= TA 2 



(A2) 



with k x the Bloch wave vector. From these boundary 
conditions we extract the relation (A = [e 2 — fc 2 ] 1 / 2 ) 



A 2 = T A 2 , with T ■■ 



w z 
z* w* 



(A3) 



For nontrivial solutions of Eq. (|A3|) , the determinant of 



A 2 — (A, .B) must be zero, i. e. 



det ( . _ ik „ ) = 0. 



(A4) 
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Evaluating the determinant gives the dispersion relation 
cosk x = ${we- iX } = cos(6 t + X)/\t\, (A5) 
with 1/w = t = \t\e iSt . 



Appendix B: Crossing points for unequal barrier 
and well widths 



Subtracting the second equation from the first one in 
(B2) gives 



2u:- u 2 (l-2W b )=ir 2 ( ± 



\ 2^'], <B3) 



wl 



Suppose a solution (e, k x — 0, k y ) of the dispersion 
relation rtTJ) is known for which the derivative de/dk y at 
a certain k v value is undefined; then this k y value can 
be a crossing point. The condition for such a solution is 
sm(\W w ) = sin(AW(,) = 0, and coa(\W w ) = cos(AW b ) = 
±1, which entails 



AWu, = jn, 

AW b = (j + 2m)n, 

with j and m integers. Explicitly we obtain 

[{e + uW b ) 2 -kl]Wl = (^)\ 



[(e - u(l - W b )f 



kl]W 2 = ((j+2m)Tr) 2 . 



(Bl) 



(B2) 



from which the corresponding value of the energy e can 
be extracted. Substituting this value in the first of Eqs. 



( B2 ) one obtains 



j 2 (j + 2mf 



2u \ Wl 



k yj m = ±[(e j<m + uW b f - (jn/W w ) 



W 2 

21 1/2 



(B4) 
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